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The problem of the quantitative degradation of the performance of a quantum computer due to 
noisy unitary gates (imperfect external control) is studied. It is shown that quite general conclusions 
on the evolution of the fidelity can be reached by using the conjecture that the set of states visited 
by a quantum algorithm can be replaced by the uniform (Haar) ensemble. These general results 
are tested numerically against quantum computer simulations of two particular periodically driven 
quantum systems. 
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I. INTRODUCTION 

In the last decade quite a few studies have been carried 
on which focus on the behaviour of an imperfect quantum 
computer (q.C.) 0. Although the sources of imperfec- 
tions are in general specific to each particular physical 
implementation, they can be grouped in two categories 
0; decoherence and unitary errors. 

Decoherence consists in an unwanted interaction be- 
tween the Q.C. and the surrounding environment, due to 
imperfect isolation |3( ■ This coupling causes the two par- 
ties to become entangled, so that the state of the Q.C. 
alone is no more described by a pure state but by a den- 
sity matrix. Its consequences on the results of a compu- 
tation, which prompted for the development of error cor- 
rection strategies, have been studied from the very early 
days of quantum computing 4] ; active corrections (error 
correction codes) are already sufficiently understood to 
be covered in standard textbooks [seeQ, chapters 8-10], 
while passive methods (decoherence free subspaces) have 
been developed more recently [sec 5, for a review]. 

Even in a Q.C. perfectly isolated from the environment 
the computation can still be affected by unitary errors; 
the state evolution remains coherent, but the algorithm 
is slightly modified. Unitary errors arise in at least two 
different contexts: they may be due to an imperfect im- 
plementation of the operations by means of which the 
computation is performed (noisy gates, N.G.) or to a 
residual Hamiltonian leading to an additional spurious 
evolution of the Q.C. memory (static errors, S.E.). 

The lack of knowledge about the parameters determin- 
ing the unitary errors implies some level of uncertainty 
for the results of the computation, that is, unitary er- 
rors can be thought of as an effective decoherence mech- 
anism. While true decoherence and N.G. are to a certain 
extent similar in their effects, S.E. differ because time cor- 



relations cannot be neglected. This is particularly clear, 
for instance, in ensemble computing, where the term in- 
coherence is used 0: even if the action of S.E. is de- 
scribed by completely positive superoperators, it can be 
corrected by locally unitary "refocussing" operations. 

The effects of imperfections can be summarised by 
means of a relation between an "imperfection intensity" 
e and the timescale t& for the degradation of some char- 
acteristic quantity linked to the computational process. 
The usual choice for this quantity is the fidelity, de- 
fined as the squared modulus of the overlap between the 
state | ip ) = U\ tp ) of the quantum memory after an ideal 
algorithm and the corresponding state = J7* | ip° ) 

after an imperfect evolution of some type, 



(1) 



There is, recently, a rising interest in the problem of 
S.E, especially in connection to different signatures for 
integrable/chaotic dynamics [1,0. The focus in this ar- 
ticle will however be on N.G, which are less explored in 
the literature. The first investigations were based on sim- 
ulations of Shor's algorithm 10]. In the seminal paper 
for the ion trap computer [111] the authors note that the 
quantum Fourier transform is quite robust with respect 
to an imperfect implementation of the laser pulses. A 
more systematic but still completely numerical analysis 
of the evolution of the fidelity in an almost identical setup 
was performed in 0. Although only Hadamard gates 
were considered to be noisy there, the empirical result 
captures the general behaviour td ~ e~ 2 . 

Grover's algorithm ^| was studied in with the 
help of a phenomenological model of probability diffu- 
sion, but the dependence of the model parameters on 
the number of qubits and e was not determined. Con- 
clusions similar to those in were reached also in [lij] . 
where again the test algorithm was Grover's one 1 , and 
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1 Indeed, the result stated by the authors in Il4l is that there is an 
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in a number of articles about the simulation of quantum 
chaotic maps 0, 0, 0, 0, 0] . These results show that 
the timescale for the degradation of the fidelity of a Q.C. 
decreases only polynomially with noise and system size 2 . 

The aim of this article is to perform a quantitative 
analysis of the effects of unbiased N.G. in a Q.C. running 
a quantum algorithm without measurements (i.e. a multi- 
qubit unitary transformation) . If G is an elementary gate 
accessible to the Q.C, its noisy implementation can be 
written as G* = £ G, where £ is a unitary error operator 
(in the following the symbol * refers to the presence of 
imperfections); so, an ideal quantum algorithm is U = 
Y\ k Gk and a noisy one is U* = Y\ k £kGk- The adjective 
"noisy" implies that there is no correlation for the error 
intensities of £k and £k'^k, an d the adjective "unbiased" 
that for an error operator £ = 1 + 0(e) the average 3 over 
different realisations gives (£) = I + 0(e 2 ). 

It will be shown that, in this "unitary model", the 
degradation of the fidelity for N.G. depends only on the 
spectrum of the errors, on the size of the Q.C. and on the 
number of gates in a given algorithm. Note that in an 
algorithm not involving measurements, the mean value of 
the fidelity is a directly measurable quantity; in fact, due 
to the previous considerations, it is sufficient to run the 
algorithm forward for N g /2 gates, and then to apply the 
adjoint algorithm to go back to the initial state (which 
is known and which can be always chosen as an element 
of the computational basis). The average fidelity for N g 
gates is then exactly the probability of measuring this 
initial state, which can be done at each fixed precision 
with a fixed number of measurements. 

The theoretical results are checked against numerical 
simulations based on two specific quantum algorithms 
calculating the evolution of a periodically driven quan- 
tum system (see appendix [XJ, the sawtooth map Q and 
the double- well map [l6j|. The former is a purely uni- 
tary transformation, while the latter uses an ancilla and 
intermediate measurements. These numerical examples 
concern thus an imperfect Q.C. used to simulate the evo- 
lution of another quantum system (with the Q.C, in turn, 
simulated by a classical computer) . Note that the choice 
of a particular set of primitive gates might affect the nu- 
merical checks but not the theoretical results, since the 
set is left unspecified there. 

The proposed model and the numerical results suggest 
that a Q.C. subject to a non-trivial computational task 



shows a sort of "universal" (algorithm-independent) be- 
haviour for the fidelity degradation. The body of the ar- 
ticle is organised as follows. Section[n]recaris the bounds 
to the fidelity degradation posed by the unitarity of er- 
rors 4 . Section II I II introduces a detailed model for this 
degradation, whose predictions are tested numerically in 
section IIVI The final part contains a discussion of open 
questions and limitations of this model. Various appen- 
dices explain the details of the calculations, of the imple- 
mentation of quantum algorithms and of their numerical 
simulations. 



II. FIDELITY DEGRADATION INDUCED BY 
UNITARY ERRORS: A SIMPLE BOUND 

As already said, a common approach in the study of 
unitary imperfections is to analyse their effects by means 
of the fidelity (equation , which can be seen as the 
probability of remaining on the state selected by the ideal 
algorithm. In a unitary error model a simple bound holds 
for its degradation; in fact, by introducing the squared 
norm of the difference of the state vectors one finds 

H1 1>) - | ^* >|| 2 = 2(1 - | ^ )) > 2(1 - v?). 

However, since / does not depend on a global phase 
change on \ip*) (i.e. [/* — ► e lx U*), while ||| -0* ) — | 
does, one can take the minimum over x- I n this case, the 
previous bound can be shown to become an equality: 
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1p* 



= 2(1-^/7). 



(2) 



The squared norm can in turn be bounded by the norm 
of the operator implementing the algorithm. If the op- 
erator is diagonalisable this norm is of course the largest 
eigenvalue (in modulus). Introducing 



|A||= max ||4 Ml 

\<t>) 



(3) 



which is normalised to 1 for unitary operators, it is im- 
mediate to see from definition [3] that 



< \\U-e lx UJ 



(4) 



The inequality \\UiU 2 - WiW 2 \\ <\Wi- Wi|| + \\U 2 - 
W 2 \\, which holds for unitary operators, implies that the 
majorisation chain can be continued with 



error threshold at e = 0(n q ~ 2 N~4 ), where n q is the number of 
qubits in the Q.C. and N = 2 n i is the number of levels. However, 
given that Grover's algorithm consists of 0(y/N) repetitions of a 
basic cycle, with 0(n q ) N.G. per cycle (only the Hadamard gates 
are perturbed), it turns out that the overall number of N.G. is 



No 



0(VNn q ), so that e 2 = 0(1/N g ), i.e. t d 



2 Although there are other quan titie s which are exponentially sen- 
sitive to the number of qubits ll :1 lid Il9l . 

3 Satisfying this condition corresponds to tuning the gate imple- 
mentation in order to eliminate the systematic part of the error. 
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\l-e lXk £ k 



e^£kGk\\ (5) 

del \ - 



4 This bounds are valid when the Q.C. evolution is completely uni- 
tary, i.e. when the executed algorithm does not use measure- 
ments before the end of the computation. 
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where \ an( i the x^'s are independent variables, since is on the small error limit, the goal is to find the value 

of a for / ~ 1 — aN g e 2 (there are however arguments 
which support the extension to an exponential law). The 
generic limit in equation [S] can be specialised for the case 
of N.G. with the value = e/8 (see formula lOT]) : 

N 2 e 2 
l-/< 9 



min x — e l ( x ~^k Xk ' > \\ = 0. In this expression, the def- 
inition of c*fc = min^ III — e**£jd| has been introduced. 
This quantity depends only on the spectrum of the er- 
ror operator £k- This is easily shown by moving back 
to definition and writing the generic vector | (ft ) over 
an eigenbasis for £k] if the corresponding eigenvalues are 
e zXkj (Xkj € R since the errors are unitary), one obtains 

<r*fe = min max (I — e lx £k)\ (ft}\\ 

X I <t>)=Y^j a kj I 4>kj ) 



= mm max 

X a kj 



mm max 

x i 



1 



i(x+Afei) 



= 2 min max 

x i 



sin 



X+^fcj 



This already shows that only the eigenvalues e lXkj of 
£k are relevant. Introducing finally the definition 



(6) 



where the average for is on the type and the intensity 
of the errors, it is possible to conclude that 
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k 



= ^N 2 . 



(7) 



Therefore, using equation [2 0] and and multiplying 
by 1 + v7 < 2, one obtains 



1 - / < ?X 



(8) 



which shows that the degradation of the fidelity in a uni- 
tary error model cannot increase more than quadratically 
in the number of gates and in the error intensity. Of 
course, this result is only a majorisation, and the true 
average behaviour can be very different. 

The quantity ^ is in general algorithm dependent, be- 
cause of the average over the type and the intensity of 
the errors, but only through the spectra of the error op- 
erators. In order to calculate a non-trivial combined 
minimisation and maximisation is required. However, if 
£k is assumed to be a small error, i.e. if all the angles 
Xkj are located in a narrow neighbourhood of radius r^, 
it can be seen by geometric means that 

^ fe = 2sin-^ < r k . (9) 



III. FIDELITY DEGRADATION DUE TO 
NOISY GATES: A PERTURB ATIVE MODEL 

The numerical simulations of the earlier articles about 
unitary errors, cited in the introduction, show that the 
degradation of the fidelity (definition ^) follows a law 
/ = exp(— aN g e 2 ), where e expresses the "intensity" of 
the imperfections (see appendix IU1 for the details about 
the error model underlying N.G.) and a is a parameter 
to be extracted with a fit. Since in this article the focus 
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As already said, the previous relation is only a majori- 
sation. In particular, equation [5] takes into account the 
worst case scenario, where all the N.G. sum up coherently. 
One could conjecture that in a non-coherent scenario for- 
mula \7\ should be modified by replacing (X^fc^fc) 2 with 

2 
9 



QZfc^*fc)> which reduces N 2 to N g ( formula IU2 



1-/ 



N g e 2 



(10) 



As already said, exception made for the numerical co- 
efficient, which is linked to the chosen error model, this 
result is completely general and can be found more sim- 
ply by noting that each noisy gate transfers a probability 
of order e 2 to the space orthogonal to \ift), and that, 
in absence of correlations, all these probabilities can be 
summed up. The law in equation I1UI for a noisy com- 
putation was first conjectured and shown numerically in 
Q, where the authors remark the fact that, although effi- 
cient non-trivial circuits in general produce entanglement 
in the Q.C. memory, it seems possible to estimate the de- 
pendence of the fidelity on e with a model where the N.G. 
errors affect each qubit independently. 

This heuristic derivation however leaves two questions 
unanswered: how to estimate the exact numerical coef- 
ficient and the fluctuations of the fidelity. In order to 
answer them, another approach will now be introduced, 
dealing with / in the limit of "small errors" . In this ap- 
proach, fk, the fidelity after k N.G.'s, is treated like a 
stochastic variable with fc-dependent distribution, with 
the constraint that fo is 1 with certainty. 

In the following it will be shown that the effects of er- 
rors can be summarised by a single quantity, the param- 
eter a 2 , which is (see also definition ID7I for cr\) the aver- 
age variance (over the error types and intensities) of the 
phases of the eigenvalues of the unitary error operators 
(the limit of small errors thus corresponds to a 2 <C 1): 



(11) 



If | ip k ) is the state of the quantum memory after k 
ideal gates, | ift k ) that after k noisy gates, and Pk = 
| ip k ) ( ift k | the projector onto the ideal subspace, start- 
ing with \tp k ) — £ k Gk \ V'* -1 ) one easily shows that 



4\^) = GfclPfc-i + Oi-Pfc-!)]!^ 1 ) 

= l^x^- 1 



I ^ l ) + 



(12) 



y/fkZil V fc > + Vl-/fc-l| Vx>, 



where e l9k is a phase such that ( ip k \ ip k ) = e l9k y/fk, and 



j ift± ) is the vector 



G fc (I-P fc _i)|^ 



fe-i ■ 



divided 
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by its norm — fk-i] note that | if)\ ) lives in the space 
orthogonal to | if> k ) . For shortness of notation let 



therefore left with 



and 



Pk 
Pk 



e -(i> k \£ k \ if) k \ 



-(^ k \£ k \^i), 



with e 1 ^ = lim at ^o(if) k \£k\ if) ). The introduction of £fc 
is motivated by the fact that in this way pk goes to 1 in 
the limit of small errors, instead of depending on a global 
phase for £k- In other words, ln£fc = - i(£fcl + Qk), where 
Qk — 0(cr*) is an Hermitian matrix. Multiplying on the 
left both sides of equation 1121 bv (if> k \£k and taking the 
squared modulus one finds 



fk — vfk-i ■ Pk + yl — fk- 



i • Pk 



(13) 



fk is therefore determined by two competing contribu- 
tions: the first term is a loss of fidelity due to the noisy 
gates moving probability out of the subspace of | ift k ) ; the 
second term represents interference from the subspace 
where | ifjj_ ) lives. 

The quadratic forms p and p, which, in absence of er- 
rors, are equal respectively to 1 and 0, have a magnitude 
which depends on ct*. By replacing \~a.Ek = *(CfcI + Qfc) it 
is easy to see that both 1 — |p| 2 and |p| 2 must be 0(er 2 ). 
Expanding the squared modulus in cauation ll3l one finds 
the following recursive relation, where the under-scripts 
stand for the orders of magnitude of the leading terms 
with respect to c* : 



fk — fk- 



i = fk-i {\pk\' + f k -i)\pkV 

0.2 



+2V7^T Vi-/fc-i»( pk Pi)- (i4) 



Now, keeping only the second order terms in formula 
1141 taking the limit of the zero order terms for ct* — ► 0, 
and summing up the partial differences, one arrives at 



N„ 



£ 

k=l 



Pk 



| 2 - 1) + 2^1- fk-i^p^ + OK). 

(15) 

The first interesting quantity to be calculated from 
equation 1151 is the fidelity averaged over many realisa- 
tions of noise, (/jv b )- Since the error intensities for dif- 
ferent £k 's are uncorrelated, the overall average splits 
into averages at fixed k. The pk term feels this aver- 
age almost only because of £k, because the vector | i\)\ ) 
is determined mainly by the previous history of the evolu- 
tion (i.e. the previous, uncorrelated errors). Since errors 
are unbiased, one gets (e~ l( ' k £k) = I + 0(er 2 ), so that 
(pk) = 0(cr 2 ) can be neglected 5 in equation IT51 One is 



5 It is not evident that this cancellation holds for s.E. or for biased 
N.G; this is the reason why 1 — (/jv g ) can be proportional to Ng 
instead of N g . The similarities between biased N.G. and S.E. have 
not been inspected so far. 



N„ 



(/^> = l + E(i^i 2 - 1 ) + °( (7 * 



(16) 



k=l 



This leading-order model returns a fidelity which is a 
function of the algorithm, because pk depends on | ip k ). 
In a different approach, one could consider the algorithm 
itself as another random variable, so that pk and pk 
would depend on three sources of "randomness" : 

the error type: the parametric form of £k, linked to the 
gate type Gk ; this allows different elementary gates 
to be affected by different types of errors; 

the error intensity: the parameter which controls the 
magnitude of the error in £k\ this allows errors for 
different realisations of the same elementary gate to 
be different; its distribution is crucial for an error 
model being unbiased or not; 

the state vectors : among which the | ip )'s depend on the 
history of the ideal algorithm and the | ipj_ ) 's on its 
noisy realisation. 

In appendix [D] it is shown how to calculate the av- 
erage value of |p| 2 when \tf) k ) is replaced by a vector of 
norm 1 randomly chosen according to the Haar (uniform) 
measure (this distribution can be generated by applying 
the unitary circular ensemble |20j to a fixed vector). It 
turns out that the distribution of |p| 2 is exponentially 
peaked around its average value; the variance is of order 
0(1/ yN), where N = T %q is the number of levels in the 
Q.C. memory and n q is the number of qubits. This means 
that, for large N, if |p| 2 is calculated on a random vector, 
the result is a value exponentially close to the mean with 
high probability: this can be rephrased by saying that 
"almost all vectors are typical" . 

This result is derived using the uniform distribution, 
but it can be argued that in order to radically change this 
picture, the distribution of the | ^>)'s should be very dif- 
ferent from the uniform one. This indeed can happen if 
the state of the system at some point factorises on the sin- 
gle qubits; given that the elementary gates are likely to be 
local to qubits, the effective ./V in this case would be very 
small (for instance, it would be 2 for 1q gates), invali- 
dating the concentration of measure plj result. However, 
it is common wisdom that efficient algorithms generate 
state sequences with non-negligible multi-qubit entangle- 
ment. So, one can be led to think that "interesting" algo- 
rithms satisfy very well the "typical vector" assumption. 
This fact could have some relation with the recent result 
that a circuit with a polynomial number of gates is able 
to approximate the statistical properties of the uniform 
ensemble, although the latter is described by an expo- 
nential number of parameters j2^ . 

For these reasons, it will be conjectured here that the 
result of equation ^| can be approximated by having 
each |pfe| 2 replaced by the average value of |p| 2 accord- 
ing to the uniform measure. This replacement, of course, 
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flushes away any dependence on the algorithm due to the 
actual sequence of the | ip ) vectors; it will be indicated 
in the following by an underline. So, the symbol ( ), like 
in a 2 = (o" 2 fe ), implies only an average over the error type 
(— ► k) and intensity. Of course, the "uniform average" 
must be taken on the actual function of p or p instead of 
calculating the function on the average of these variates; 
in other words, it is consistent, for instance, to replace p 
with zero but retain |p|. 

Due to the aforementioned assumption \pk\ 2 can be 
replaced by \pk\ 2 — 1 — Aa\ k + 0(a\ k ) (see formula lD8l 
where A —n/(i+n) is a vanishing dependence on the di- 
mension of the state space very close to 1); formula [TBI 
becomes then 



0{*l 



fc=i 



Xk I 



AN g ol- 



(17) 



natural approximation is to keep only the term scaling as 
NgE 2 in formula IT51 So, for N g » 1 formula IT51 becomes 

a 2 (f Ng ) ~ 2N 2 g a 2 -(a 2 {Kp k )) = 2N 2 cr 2 «5R 2 p fc ) - <5Rp fc } 2 ) 

As before, the average over the error realisations for 
each p k is null, therefore (SRpfc) — > 0. Then, expanding 
the real part the previous formula becomes 



a 2 (f Ng )^^-N 2 a 2 (p 2 



~*2 
Pk 



2\p k \ 



The uniform average prescription for p k consists in 
choosing | ip k ) uniformly random and | ) uniformly 
random in the orthogonal subspace. As shown in ap- 
pendix [D] the term p 2 , averages to zero, and the term 
\pk\ 2 can be replaced by \pk\ 2 — N^T a xk (formula [D9j). 
Approximating N ¥_ 1 *— > jj just like A i— ► 1, it is possible 
to conclude that 



this result agrees with formula 1 101 but for the numerical 
factor (which has however the correct order of magnitude 
since §|<cr^<f|, according to formula l?T|) and the slight 
dependence on N given by A (hard to prove numerically) . 

The second interesting quantity to be calculated is the 
variance a 2 (fw g ) ■ It is difficult to extract the exact vari- 
ance from the model in equation 1151 however, since it 
is used only to decide the number of retries for numer- 
ical simulations, an approximate solution is acceptable. 
If cr{fN g ) < 1 - (fN g ) then fk-i in formula [TBI can be 
replaced by its average value given by equation El The 
simplified model becomes then 



E[(i 

fc=i 



pk 



\ 2 -l) + 2a^A{k-l)^p k . (18) 



Forgetting for a moment the average on the algorithm, 
the variance of the fidelity can be calculated using the 
identity a 2 {y) = J2k a2 ( x k) + J2 k jth c °v(^fe, %h), valid 
for a variate y — J2k Xk - When this identity is applied 
to formula 1181 only the correlations between the terms 
with the same gate index k survive; therefore, setting 
the factor A to 1 (which it is exponentially close to), and 
confusing N g ± 1 with N g where not sensitive, one finds 



I> 2 ( 



Pk\ 2 )+^ 2 j2 k (k-iy^Pk) 
k) 



+4cr*^ fc V'fc - 1 • cov(|p fc | 2 
N g (a 2 (\p k \ 2 )) +2N? j( t 2 ■ (a 2 (Kp k )} 

+ |A 9 3 / 2 ^.(cov(|p fc | 2 ,5Rp fc )> 



Note that the three terms in the previous expression 
are proportional respectively to N g 1 e 2 , e 2 and N g 1 ^ 2 e 2 , 
where e = N g a 2 ; therefore, for a large number of gates it 
is the second term (that with the variance of #lp k ) which 
dominates. Since the parameter which expresses the dif- 
ficulty of extracting er* from a single numerical simula- 
tion is the ratio between the fidelity standard deviation 
cr(/jv ) and the average fidelity decrease 1 — (/jv ) oc e, a 



° 2 (fN g ) 



N 2 a 2 



\Pk\ 



= M N ^ 2 ) 2 = M 1 -(fN g )) 



(20) 



The ratio between the fidelity standard deviation <j(fN g ) 
and the average fidelity decrease 1— (/jv fl ) depends there- 
fore only on n q (remember that the number of gates is 
N g while the size of the Q.C. memory is N = 2 nq ). 

Two general remarks should be made on this detailed 
fidelity model. First, strangely enough, in eauation ll5l it 
is the first term which determines the average fidelity, but 
the second one which governs its fluctuations. If the con- 
tribution of p in relationHTjlwas null, then <J 2 (fN g ) oc N g 
and the ratio in formula QUI would drop to zero as N^ 1 , 
i.e. the fidelity evolution would be self-averaging; more- 
over, /tv s would decrease monotonically as a function of 
N g ; both these consequences are contradicted by numer- 
ical simulations (see, for instance, figure EJ. 

Second, in the derivation of the statistical properties of 
/jv s the number of gates N g was regarded as a constant. 
This does not mean that the fidelity values at different 
stages of the algorithm are not correlated. Indeed, this 
correlation can be large, because of the limit on the fi- 
delity decay discussed in section ITT1 

In the specific error model studied in this article (see 
appendix O) . a 2 depends only on the number of qubits 
which the error operators act on (1q or 2q errors). For 
R v (0 the eigenvalues are {e _? / 2 , e^/ 2 }, and one gets 
a\ = £ 2 /4 + 0(e 4 ). On the other hand, for controlled 
phase shift errors, the eigenvalues are {1,1,1, e 1 *} and 
the variance is a 2 — 3£ 2 /16 + 0(e 4 ). Thus, in general, 
(19) CT 2 _ (j^2 _|_ 0(e 4 ), with C depending only on the gate 



being 1q or 2q and £ only on the error parameter distri- 
bution. The mean value (c 2 k ) implies both an average 
on the error type and intensity, giving 



v 2 = (*i k ) = (c)(e) = 



n\C\ + n 2 Q 



12 



V, 



(21) 



where < n 1 (n 2 ) < 1 is the average number of lQ (2q) 
gates during the computation (so that n\-\-n<i = 1). Since 
d = \ and C 2 = j|, this implies 48 < (e 2 /cr 2 ) < 64. 
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Figure 1: (Color online) 

(a) Ratio {z\f \ > where / is the unaveraged fidelity and (/) 
is its theoretical value given by formula IT7I for 100 numerical 
simulations on a sawtooth map with n q = 12 and e = 0.002. 
As in the other numerical simulations in this article involving 
the sawtooth map, the number of cells is L — 2, the classical 
parameter is K = 0.04, the fidelity is observed in position 
representation and the initial state is H <g> I ® H® {nq ~ 2) \ 0); 
see appendix for further details. The same behaviour (a 
spread depending only on the number of qubits for N g S> 1) 
was found for all the numerical simulations used in figure^ - 
The largest fluctuations tend to be located below the average. 

(b) Values of Y~TJ)^ each point corresponding to 676 nu- 
merical simulations at fixed error intensity e and number of 
qubits n q ((/) is the numerical average of /). For each n q 
in {3, ... , 14}, the error intensities e = 0.002, 0.005 and 0.01 
were tried. Formula l2Q| predicts that these points depend only 
on n q . The exact dependence y/l/N, thick line, is however 
poorly followed. The thin line shows a fit with an exponential 
function aN~ b , giving a ~ 1 and b ~ 1/3. 



IV. NUMERICAL CHECKS 



The described fidelity decay model has been tested 
numerically on two algorithms simulating periodically 
driven quantum system: the sawtooth map and the 
double- well map [Tg]- These systems are usually stud- 
ied in the quantum chaos field, because in spite of their 
relative simplicity they show a large variety of physical 
phenomena, from dynamical localisation to quantum er- 
godicity, first studied on the kicked rotator model |2^| . If 
f] and 6 are momentum and coordinate operators, then 
the evolution corresponding to one map step for a map 
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Figure 2: (Color online) Decay of the fidelity due to N.G. 
during the simulation of the sawtooth map, with n q = 10 and 
e = 0.003. The simulation of one map step for such n q requires 
n g — 200 gates, with m = 45 lQ gates and = 155 2q gates. 
The value of F from formula 1171 is therefore ~ e 2 n s /59.6, 
shown by the straight line. The red (dark) curve corresponds 
to an average over 500 numerical simulations (it is hardly 
distinguishable from the straight line in this picture), while 
the orange (light) curve to a single try. See figure 0a for the 
map parameters. 



with potential V(9) is given by 



Up 



exp 




(see formula IA3I and other physics and implementation 
details in appendix [XJ . The behaviour of the classical 
counterparts of these chaotic maps is controlled by a sin- 
gle parameter K — kT, which in the following simula- 
tions is set to 0.04. For this value, the double-well map 
presents a mixed phase space, with a large island of in- 
tegrable motion in a chaotic sea; the sawtooth map is 
chaotic for every K > 0, but the cantori regime extends 
up to K ~ 1, so that the presented numerical data best 
describe a perturbative regime. In both cases, it is not 
evident that the success of the uniform measure approx- 
imation depends on "chaos" . 

The algorithms corresponding to the simulation of 
these maps contain the Fourier transform as a basic in- 
gredient for changing the representation from coordinate 
to momentum and back. Since the simulated maps have 
periodic impulsive kicks, it is possible to write the overall 
evolution as a piecewise diagonal unitary transformation 
in the appropriate representation. The set of elemen- 
tary gates which are used for the algorithms implemen- 
tation contains single and double qubit phase shifts and 
Hadamard gates (see appendix[U]for more details). Since 
these algorithms are periodic, the results are usually ex- 
pressed in terms of the number of gates n g per maps step. 
In the following the variable t will be used when refer- 
ring to the "map time" (the number of map steps). It 
is immediate to specialise the previous results for a map 
simply by replacing N g with n g t. The effective decoher- 
ence parameter T is defined by / ~ 1 — Ft. 
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As a first check, the theoretical predictions were com- 
pared with the results of the simulation of a sawtooth 
map. This map is similar to the double- well map but it 
is free from the complications of using auxiliary qubits. 
Figure is a check of formula [3UJ it shows that indeed 
the ratio <j(/jv )/(l — (/jv b )) depends only on n q and is 
exponentially small with it, but the theory predicts a de- 
crease like N~ x / 2 while iV" 1 / 3 is observed. A possible 
explanation could be that the vectors \ ip\) do not ex- 
plore the complete (N — l)-dimensional space orthogonal 
to | i\) k ), but only ~ N 2 / 3 dimensions. 

The previous ratio characterises the difficulty of fit- 
ting r from a single numerical simulation, and can be 
used to estimate the number n S i m of tries necessary to 
reduce the statistical uncertainty, for instance setting 
^== \-(f) ~ This prescription implies a number of 
simulations n s j m ~ 10 4 2 - ™'; however, it is better to use 
the numerical values from figure 0-b, because the theo- 
retical formula gives an underestimate. In the remaining 
of the article n S i m varies between 50 to 500. 

An example of the degradation of the fidelity for the 
sawtooth map is given in figure [21 where it can be clearly 
seen that a single numerical simulation can give a fidelity 
which is not monotonic and correlated over many map 
steps. The prediction T — An g a 2 , given by formula 1171 
is compared with numerical simulations, again for the 
sawtooth map, in figure|31 It turns out that the theoreti- 
cal value is an upper bound for the the numerical fidelity 
decay constant, and that, not surprisingly, the agreement 
is better for a large number of qubits and small errors; 
for n q > 10 1 and e < 10~ 2 , the discrepancy is within 5%. 

The simulation of the fidelity decay on the double well 
map presents a new effect: beyond "standard" fidelity 
fluctuations, there are some rare "anomalous" fluctua- 
tions, like that in figure with the fidelity dropping 
abruptly to almost zero. This is due to the fact that the 
algorithm implementation, for n q > 4, uses an auxiliary 
qubit (an ancilla), which is reinitialised to | ) every time 
it is reused; the reinitialisation implies a measurement, 
which, when the evolution is affected by noise, can select 
the "wrong" result (i.e. | 1 )), thus in practice completely 
destroying the computation. These catastrophic events 
become more and more frequent while e increases. 

However, when the measurement of the ancilla gives 
the "good" result, it turns out that the degradation of 
the fidelity is slowed down. This cancellation of the (un- 
wanted) evolution is known as Zeno or watchdog effect 
|24| ; recently it has been shown that this inhibition of de- 
coherence can arise in more general contexts, where the 
essential ingredient is a strong coupling to the environ- 
ment and not its trivial dynamics 25j . In the simulation 
of the double- well map it is indeed the ancilla which plays 
the role of the environment. 

The repeated observation of the ancilla can be thought 
of as an error correction strategy in the Zeno regime [26| , 
because the Q.C. memory state is partially rectified even 
when the measurement gives the expected result. The 
Zeno effect is known to be linked to error correction from 
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Figure 3: (Color online) 

(a) Ratio F/F t h versus number of qubits n q for the sawtooth 
map. r is the (numerical) fidelity decay constant and T t h is 
its theoretical average, given by formula |l7l Each point cor- 
responds to an average at fixed n q and e with 500 numerical 
simulations. Yellow (light) points correspond to large errors 
(e ~ 10 _1 ) while red (dark) points to small errors (e ~ 10~ 3 ). 
The solid line is an eye-guide for e = 0.004. See figure0a for 
the map parameters. 

(b) Same as picture 0a, but the independent variable is the 
error intensity e, while the yellow to red (light to dark) transi- 
tion corresponds to an increasing number of qubits. The solid 
line is an eye-guide for n q — 11. See figure 0a for the map 
parameters. 

(c) Same as picture 0-a/b, but the dependence on the error 
intensity e and the number of qubits n q is shown at the same 
time. Each coloured box corresponds to an average at fixed 
n q and e with 500 numerical simulations. The yellow to red 
(light to dark) transition corresponds to an increasing value 
of the ratio, which remains however in general limited by 1. 
Missing boxes are due either to combinations taking too long 
to simulate or to fidelities dropping too fast. See figure 0a 
for the map parameters. 
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Figure 4: (Color online) Comparison of a typical and anoma- 
lous decay of the fidelity for the double- well map, with n q — 5 
(plus one ancilla) and e = 0.01. Both curves correspond to 
a single numerical simulation. The anomalous drop is due to 
a measurement of the ancilla giving the "wrong" result (note 
the change in the vertical scale). As in the other numerical 
simulations, unless otherwise stated, the number of cells is 
L = 2, the classical parameter is K = 0.04, the centre of the 
potential well is a = 1.6, the fidelity is observed in position 
representation and the initial state is H ® I ® H®(" q ~ 2) \ 0) 
(see appendix lAt. 



the very early days of quantum computation |27| ; it has 
been suggested as a stabilisation strategy by Shor [T^j m 
his famous paper on prime factoring, in a context identi- 
cal to the current one (i.e. syndrome measurement with- 
out a recovery circuit). 

Of course, a direct observation of a fidelity jump, like 
that in figure is not possible on a real Q.C. It is true 
that the result of each ancilla measurement is accessible, 
so one could decide to purge a set of real experiments of 
those instances which showed a "faulty" reinitialisation; 
another option could be not to reinitialise the ancilla at 
all. In any case, the only meaningful quantity is the 
average fidelity. 

The numerical simulations presented in this article are 
not purged of the unwanted reinitialisations. This slow- 
down is illustrated in figure where the numerically 
found value of T is approximately 10% to 20% smaller 
than the prediction of formula 1 171 with a preference for 
~ 15%. Note that for n q — 3 the algorithm implementa- 
tion does not use ancillae, and the result is in agreement 
with what was found for the sawtooth map (figure |3J • 
"Wrong measurement" events, like that shown in figure 
0] become more and more frequent in the averaging en- 
semble for large error intensities. 



V. CONCLUSIONS 

In this article the degradation of the performance of a 
Q.C. affected by N.G. is investigated, both analytically 
and numerically, by studying the fidelity of the noisy 
states with respect to the ideal state (definition [JJ : the 
average value of the fidelity after a fixed number of el- 
ementary gates characterises the "precision" of the Q.C. 
and can be efficiently extracted from a set of experiments 
via forward-backward evolution. The numerical tests are 
performed using two algorithms that simulate the saw- 
tooth and the double- well maps. In the following n q is 
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Figure 5: (Color online) 

(a) Ratio Y/Tth versus number of qubits n q for the double- 
well map. r is the (numerical) fidelity decay constant and Y t h 
is its theoretical average, given by formula [T7l Each point cor- 
responds to an average at fixed n q and e with 200 numerical 
simulations; at least 20 map iterations can be simulated be- 
fore / drops below 0.9. Yellow (light) points correspond to 
large errors (e ~ 10 _1 ) while red (dark) points to small errors 
(e ~ 10 -3 ). For n q = 3 the algorithm does not use ancillae. 

(b) Same as picture Qjja, but the independent variable is the 
error intensity e, while the yellow to red (light to dark) tran- 
sition corresponds to an increasing number of qubits. Most of 
the points for n q > 4 (filled circles) lie in the 80% — 90% re- 
gion, with a preference for ~ 85%. Empty circles correspond 
to the case n q = 3 (no ancilla). 

(c) Same as picture |SJ-a/b, but the dependence on the error 
intensity e and the number of qubits n q is shown at the same 
time. Each coloured box corresponds to an average at fixed 
n q and e with 200 numerical simulations. The orange and red 
(dark) boxes correspond to a T which is ~ 15% smaller than 
its theoretical value due to the reinitialisations of the ancilla 
(Zeno effect). For n q — 3 the algorithm implementation does 
not use ancillae. Missing boxes are due either to combinations 
taking too long to simulate or to fidelities dropping too fast. 
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the number of qubits in the Q.C. and N = 2"« is the 
number of states. 

The degradation of the fidelity is shown to be deter- 
mined basically by of, the mean variance of the phases 
of the eigenvalues of the error operators (definition ; 
stated in this way, this result is valid for every choice of 
elementary gates. In a realistic case, where the Q.C. is 
endowed with a particular set of elementary operations, 
each characterised by a specific error operator distribu- 
tion, this quantity is algorithm-dependent only through 
the relative frequencies of the elementary gates; so, of re- 
flects the effective decoherence inherent to the Q.C, and 
a careful analysis of the results of different algorithms 
would allow in principle to determine the proper for 
each gate type. 

The noise-averaged fidelity after N g elementary gates 
is found to be (f Ng ) = 1 - j^Ng^ + 0(af). The flue- 
tuations of the value of the fidelity after a fixed number 
of gates are such that o(/)/(l — (/)) depends only on 
n q and is exponentially small with it. It is also found 
that the reinitialisation of the auxiliary qubit during the 
computation for the double-well map reduces the fidelity 
decay rate of ~ 15% (this is known in literature as an 
error correction in the Zeno regime). 

Note that a particular case of distribution of the error 
operators is the one which does not depend on the gate 
type (this case bears the additional bonus that there is 
no dependence at all on the algorithm). In this scenario, 
an ensemble of unitary errors is equivalent to the more 
commonly studied case of a non-unitary error channel, 
like for instance the depolarisation channel. 

The rather good agreement between the quantitative 
predictions of the theory and the numerical results for the 
double-well and sawtooth maps suggests a general valid- 
ity for the assumption on the distribution of the states 
visited by an algorithm. It is clear however that there 
exist trivial algorithms which do not follow the predicted 
behaviour for some choice of their initial state. The ques- 
tion is then to characterise the class of algorithms which 
satisfy well enough the predictions, and, in particular, 
to understand whether this is linked or not to the fact 
that the algorithms used as examples exploit the Q.C. in 
order to simulate another quantum system, showing er- 
godic properties. Due to the robustness of the prediction 
for of for state distributions close to the uniform distri- 
bution, the conjecture is that the class is not limited to 
this type of algorithms. 

It would be interesting to extend the ideas presented 
in this article in order to include measurements of aux- 
iliary qubits; as shown by numerical simulations, mea- 
surements of qubits supposed to be, in absence of errors, 
in a known state can slow down the fidelity decrease. In 
perspective, this could allow to deal with error correction 
codes, a subject which has not been treated here. 

The numerical simulations in this article were per- 
formed using a freely available implementation of the 
quantum programming language described in |28j . The 
author acknowledges useful discussions with A. Pomer- 
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Appendix A: HAMILTONIAN MAPS AS 
QUANTUM ALGORITHMS 

The goal of this appendix is to show the general ap- 
proach to build a quantum algorithm for a quantum 
Hamiltonian map derived from a classical unidimensional 
kicked map. A classical kicked map is the stroboscopic 
observation of the phase space of a system which is af- 
fected by an impulsive conservative force given by a pe- 
riodic potential at regular intervals in time and which 
evolves freely in the meanwhile. If the observation oc- 
curs just before the kick, the map connecting the points 
in the dimensionless (9, rf) phase space can be written, 
without loss of generality, as 

f v = v - fc w / Ai \ 

\6 = 8 + Tfj (mod 2tt) { ' 

In this expression, T represents the (dimensionless) 
time interval between two kicks, 2ir is the period of the 
potential V{9) and A; is a parameter which governs the 
intensity of the potential. It is easy to show that, if the 
rescaled variable r\' = Trj is used instead of 77, the map 
depends on one parameter only, that is K = kT. In the 
map lAll the 9 range is finite, while the momentum 77 can 
vary arbitrarily. However, the transformation 

(6>, 77) -> (6 1 i 1 + 2irL/T), L e Z, (A2) 

is a symmetry for the system; this invariance identifies a 
natural cell, with an extension of 2tt in both directions of 
the phase space, and L can therefore be associated to the 
number of cells. The time-dependent classical Hamilto- 
nian driving the evolution described by the map I All is 

2 +°° 
H(9, V ,t) = ^- + kV(9) J2 S(t-jT). 

j=—oo 

The map quantisation can be accomplished by replac- 
ing the dynamical variables 9 and r\ with the correspond- 
ing operators 6 and 77. In the 9 representation, the mo- 
mentum operator becomes, as usual, 77 = —i-§g- The 
quantum equivalent of the map IA1I is the evolution op- 
erator calculated for a time interval T (which includes 
always one kick), known as the Floquet operator (Uf)- 
Given that the action of the potential is impulsive, the 
integration of the quantum Hamiltonian is very easy: 
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Therefore, in the quantum case the dynamics is deter- 
mined not only by the product K = kT, but also by the 
"magnitude" of h (for instance, by the ratio k/H). 

The simulation of the quantum map on a finite Q.C. 
implies some type of discretisation, which introduces an 
additional parameter, namely the number of available 
levels, N = 2 n <>. The easiest choice is to associate the 
computational basis states | i ) to some eigenvalue 0i . In 
this way the "position" operator 9 becomes diagonal in 
the computational basis. In the algorithm described in 
this article the following choice holds: 



jj^ + Se). 



(A4) 



In this way the position eigenvalues are equidistant; 
Se is an offset which fixes the exact 2ir interval where 
the eigenvalues are located. The choice Sg = (1 — N)/2 
corresponds to — 7r < 9i < tt, with perfect symmetry 
around 9 = 0. The periodicity condition 9 = 9 + 2tt 
implies that the eigenvalues of the momentum must be 
quantised in units of h. The finitcness of the Q.C. memory 
allows only for a finite number N of eigenvalues of the 
momentum to be simulated. It is therefore natural to 
resort to a cyclic condition, that is assuming that the 
range of the eigenvalues corresponds to an integer number 
of cells (see the transformation |A"2|) : 



Ar] = Nh = 2irL/T 



2ttL 



(A5) 



The last step in the construction of the algorithm is the 
fact that the quantum Fourier transform T "exchanges" 
the momentum representation with the position repre- 
sentation. This means that if, at some time during the 
simulation, the Q.C. memory's state is \ip) = X)j a jli) 
and it is associated to the simulated system's state 
Ylj a j where the | 9j ) 's are the eigenstates of 9 with 
eigenvalue 9j , then | -0 ) = J2j bj I j ) must be associ- 
ated to the simulated system's state bj \rjj), where the 
| rij )'s are the eigenstates of t) with eigenvalue r)j = jh. 
In other words, the Floquet operator Up (see lA3|) can be 
simulated by a quantum circuit Q whose action is: 



Qv\f, 
Qe\j\ 



= exp 



h 

Jim- 



N m 



(A6) 
(A7) 

(A8) 
(A9) 



It is now useful to define an operator U (/3,p) as 

U(p,p)\j) =exp(2ni/3f)\j). (A10) 

The circuit Q v a cting in the momentum representation 
(see equation IA7I can then be simplified by using the 



definition IA10I and replacing the eigenvalue rjj with jh 
and the product Th with the help of expression I A5I 



J 



= U 



-L 

2N' 



\3 , 



(AH) 

In this expression the parameters T and h have been 
replaced by the number of levels N and the number of 
cells L: the free evolution is thus completely geometrical. 
The (dimensionless) quantum action h is proportional to 
L/N, therefore the semi-classical limit can be attained 
by taking L/N — > 0. 

The circuit Qg (see equation IA8|) is easy to manage 
when the potential V(9) of the map I All is a polynomial 
in 9 (let p be its degree). One can in fact introduce the 
polynomial V{x + Sg) = YZn=v fm,x m and set 



W> = i-(77(j 



Se), 



where 9j has been replaced with expression IA4I The 
circuit Qg (see equation IA8|) can then be broken down 
into circuits of type U (/3, m) (see the definitions I A 1 01 and 
IA5I and remember that K = kT): 



Qe\j] 



exp 



% 2ttL/N J lJ) 



II" 



-Kf m (2n) 
LNP- 



p-2 



(A12) 



Note that the term m = 0, which gives only a g lobal 
phase, can be neglected. For the double-well [l6j map 
(p = 4) the polynomial is: 



Hv) 



where a is an additional parameter which locates the cen- 
tres of the two wells. The coefficients f m are therefore 

U = l, f 3 = iS e , f 2 = 6S 2 e -2a 2 , h=A5e(5 2 e -a 2 ). 

For the sawtooth |9] map (p = 2) the polynomial and 
the coefficients are respectively 



V{y) 



l 



Summarising, it has been shown that the circuit Q 
(equation IA6|) can be written as a composition of cir- 
cuits for the Fourier transform F (equation IA9(I and for 
the operators U (/?, to) (equation IA10|) . The parameters 
which govern the simulation are the classical parameter 
K (integrability — > chaos) , the number of simulated cells 
L and the number of available levels N (L/N — > be- 
ing the semi-classical limit) in addition, of course, to the 
functional form of the potential V(9). 

The state of the Q.C. memory after each application of 
the circuit Q corresponds to the coordinate representa- 
tion for the simulated system, but trivial changes allow 
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the simulation in the momentum representation. It is 
well known that the circuit for T can be implemented effi- 
ciently [29j . Appendix IbI shows that also U (0,m) can be 
implemented efficiently; the conclusion is therefore that 
the circuit for the simulation of a quantum kicked Hamil- 
tonian with a polynomial potential can be implemented 
efficiently on a Q.C. 

There exist, of course, interesting potentials which are 
not polynomial; for instance, the kicked rotator map has 
V(9) = cos6>. For this case, two approaches are known. 
In |30( an auxiliary register and a number of additional 
ancillae is used for calculating | 6 ) | ) — * | 9 ) | cos 9 ) with 
a finite number of digits in the mantissa; the auxiliary 
register is subsequently used to implement Qg. If the 
size of the mantissa for the calculation of the cosine is 
0(n q ), the circuit depth is 0(n^) and the number of an- 



cillae is 0(n q ). In |31| another approximated method, 
with circuit depth 0(kn q ), is introduced. This method 
is very interesting because no auxiliary qubit is required; 
however, the precision with which Qg is implemented is 
0(k), so the approach is optimal only for small k values. 



Appendix B: A QUANTUM CIRCUIT FOR 
EXPONENTIATION 



The goal of this appendix is to show a procedure for 
building an efficient quantum circuit implementing the 
unitary transformation 



U(0,p)\x\ 



(Bl) 



In expression IB II the exponent p is a positive integer 
while f3 is a real coefficient. | x ) is the x-th element of the 
computational basis, therefore IA (f3,p) is diagonal in this 
basis. The circuit is implicitly parametrised by the num- 
ber n q of qubits in the register which hi p) operates 
on. The first step is to write the integer x labelling | x ) 
as a binary string, x = X^=o ■ By replacing this 
expression into the phase of definition IB1I it is possible 
to rewrite it as the product: 



^2rrif3x p _ 



n 

ji ■ ■ -j P 



27ri/3aj 1 ...a 3p 2 J i + -- +3 P 



Since the coefficients aj are binary digits, the product 
a,j 1 . . . a j is zero (i.e. the phase is trivial) unless all the 
(Zj's are equal to 1. Therefore each factor in the previous 
expression corresponds to a multi-controlled phase gate 
in the circuit for U ((3,p) of definition IB II This gate acts 
on the qubits selected by the set W(ji . . . j p ), which con- 
tains the values of the indexes ji, ... j p with neither rep- 
etitions nor order. Introducing the notation Cwisp) for a 
multi-controlled phase gate applying the phase exp(27ri</7) 
to the qubits in the set W, one obtains 



u(j3, P )= n c 



W(j x ...j p ) 



092 



3i H Hp) 



(B2) 



Up to now, this method follows exactly the procedure 
described in [16j . The set W does not contain duplicated 
elements, so that its cardinality is smaller or equal to p 
(but it is always positive). Since neither the sum j\ + 
■ ■ ■ + j P nor the set W(ji . . . j p ) depend on the order of 
the indexes, all the gates concerning the same indexes 
and differing only in their order can be collected into a 
single gate. Let 



P = {JC{0,. 



1}|J^0A|J| <p}. 



(B3) 



where . 



1 are the possible values for the qubit 



indexes. Then expression IB2I can be rewritten as 

«om = n n c w{n... h) {^ h+ '" +jp ) 



JeP ji-jp 

W(ji...j p )=J 

JeP h—jp 

WUl-3 P ) = J 



(B4) 



where in the last line the gates with the same J have been 
compressed into a single gate. Note that two gates act 
on the same qubits if and only if their W sets are equal. 
To proceed, one needs to define the set Gk{n) of the 
partitions of n objects in exactly k non-empty subsets. 
A partition g 6 Gk{n) is then a sequence of k positive 
integers {gi ■ ■ ■ gk} such that their sum is equal to n. 
Given a set J and a partition g G G\j\ (p), the pair (J, g) 
corresponds to a subdivision of the set of p-tuples whose 
W set is J. It is then possible to apply the replacement 



E 



W(ji...j P )=J 



E E 

seGijlfa) 3i— ip 

W(j!. ..]„) = J 

31 H Hp=E 9iJi 



(B5) 



For each pair (J, g), the sum ji + ■ ■ ■ + j p is fixed, 
therefore the last summation on the right of expression 
B5I can be replaced by its multinomial weight 



gi'----g\j\< 

Putting all together, one finally obtains that the operator 
Li {f3,p) can be replaced by a product of multi-controlled 
gates Cj(ip), with a bijective correspondence between the 
gates and the sets JeP, where 



<P = P 



E 

s6G|j|(p) 



pi 



. 2E 9iJi 



9i- ■ --9\J\ 



(B6) 



This gate collection method presents the obvious ad- 
vantage that the circuit depth decreases. In fact, if one 
gate is built for each possible index combination ji . . . j p , 
the total number of gates is trivially (n q ) p . If, on the 
other hand, the previously described compression is ap- 
plied, one obtains as many gates as the cardinality of the 
set of partitions P (definition IB3(1 . i.e. 



min(n g ,p) 

E 

fc=i 



(B7) 
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For n q < p this sums to 2 7lq — 1. However, it is more 
interesting to consider the limit n q p; approximating 
the sum with the help of Stirling's formula one obtains 



IPI 



It is easy to check that this approximation is already 
good for n q — 5 and p — 4, where the savings 1/p! corre- 
spond to ~ 95% of the gates. It should be kept in mind 
that the calculation of tp in expression [HU can be numeri- 
cally critical, because for large n q the integer multiplying 
(3 grows exponentially as 2 prlq . Since this phase term is 
to be taken modulo 1, the least significant bits of the 
representation of are to be handled carefully 6 . 

In practice however, when U (0,p) is used as a block of 
a larger circuit, like lAlll or IA12I (3 is also exponentially 
small. For p > 2 the worst case comes from circuit IA11I 
where (3 = 0(2-™"), so that the phase tp is 0(2 n °<*-V). 
If j3 is stored as a floating point number, one loses p — 
1 bit of the mantissa of tp for every qubit added. For 
a double precision floating point 7 this problems arises 
around n cr jt = 18 for the double- well map and around 
"crit = 53 for the sawtooth map. This is beyond or at the 
limit of the possibilities of current Q.C. simulators. In any 
case, the study of the real degradation of a computation 
due to this loss of precision in the most critical gates for 
n q > n cr ; t is not trivial and has not been performed yet. 



Appendix C: ERROR MODELS AND BOUNDS 
TO THE FIDELITY DEGRADATION 



In all the numerical simulations of this article, it is 
supposed that the Q.C. is endowed with the following set 
of hardware operations: 



' l l * 
.1-1, 



R<j> — ( 1 e >*) CR<f) — 



where <f> £ K is a phase (indeed, it can be only a multiple 
of some "atomic" phase, like 27r/2 32 , which mimics the 
limited precision of the classical control system) . H is the 
Hadamard transformation, while R^ and CR$ are phase 
shifts and controlled phase shifts respectively. In general, 
the controlled operation CA corresponds to the matrix 
( 1 A ) . It is sometimes useful to express the elementary 
gates as rotations i?^(a) of an angle a around an axis r): 



Rfj(a) 



cos 



(!H sin (f) 



r\ ■ a . 



With this notation, setting xz = [x + z)jy2, one 
finds H— e l % R^ z (tt) = R^ z (n) and R$ = R z (4>) = 
Rz{4>), where '=' means equivalence modulo a global 
phase (which can be neglected). Note, however, that 

CRj, = e l i[R z {^)®l]oCR z {(j)) jk CR z {<j)), so that CR^ 
and CRz(4>) are inequivalent even modulo a global phase. 

In the N.G. error model (see for instance 0]), each 
gate is replaced by a still unitary transformation, close in 
norm to the original gate, parametrised by a single error 
variate £ with flat probability in [— e/2, e/2]. The quan- 
tity e is called the "error intensity" , and it summarises 
the amount of noise which affects the Q.C. It is useful to 
express a perturbed elementary gate as a composition of 
the unperturbed gate followed by an error operator close 
to the identity. The imperfections for ^-rotations and 
controlled z-rotations are implemented by shifting their 
phase angle by an amount £; it is easy to show that 

noisy R^ — ► R^R^ = R^R^, and 



noisy CR^ 



CR(CR$ = 



R 4 



CR z {0 CR 4 



Each Hadamard gate is transformed into a rotation of 
an angle 7r around an axis ff randomly tilted around the 
unperturbed direction i) = xz of an angle 8 £/2. The 
perturbed gates can be generated by taking a random 
direction fi in the plane orthogonal to xz and setting 
ff = fj cos I + /t x f) sin | . Note that the signs of [i and £ 
are not independent. A little algebra shows that 

noisy H — ► Rp,(QH. 

Thus, the error operators £k in the N.G. error model 
are i?z(£), Rp.{£) and Ci?j, the phases of the eigenval- 
ues for 1q errors are {±|}, and those for 2q errors are 
{O3, £}. All these eigenvalues are of course specified mod- 
ulo a global phase shift. It is easy to see that this model 
is unbiased. Another approach to N.G. [32] consists in 
diagonalising the gates and perturbing the eigenvalues of 
an amount £; given that the average properties of the in- 
duced effective decoherence are determined only by the 
error operators' spectrum, it is not surprising that this al- 
ternative approach yields in general very similar results. 
The calculation of the radius of the neighbourhood of 
the eigenvalue phases (see equation |?5J is very simple; 
is always |£|/2, so that (see equation|fjJ 



The incoherent replacement (^2 k ?* 



(CI) 

(Efc<4) in 



formula [7| gives the result necessary for formula ITUI 

'C 2 \ N r 



1 - / < N g 



e L 4 s 48 y J 



6 This problem of course does not originate from the collection 
method. Factors of order 2 pn i are present also in the "plain" 
approach. 

7 Take for instance the 64 bit IEEE floating point format, which 
has 53 binary digits available for the mantissa. 



Contrary to previous works Il5l . the tilting angle for H is £/2 
and not £. This uniforms the effects of lQ gates, because all lQ 
error operators are rotations of an angle £ around some axis. 
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Appendix D: UNIFORM AVERAGES AND 
CONCENTRATION OF MEASURE 



This appendix details the evaluation of an integral like 



1(R 



(i) 



,i? w ) 



duty) Y[(i>\R {k) \i>), (Dl) 



k=i 



where the R^> are operators over an iV-dimensional 
Hilbert space, and dfi(ip) is the uniform measure over 
its unity vectors (that is, the unique invariant measure 
under the action of all unitary matrices) . If the integrand 
is expanded over some basis, then 

i(i?«...,i?M) = n<i- jdm^lWi* > 

1„ J 1- 



where the linearity of the expectation value and the fact 
that the average is over the | ip ) 's and not over the R w 's 
were used (summation over repeated indexes is under- 
stood). The measure d^(tp) over the Hilbert space can 
be turned into a measure d/J,(U) over the set of unitary 
transformations of this space. Since dfj,(tp) is by defini- 
tion invariant when U\ip) is substituted for in for- 
mula lDll with U a generic unitary matrix, d^(U) remain 
completely determined (Haar or uniform measure, corre- 
sponding to the unitary circular ensemble |20|). 

This invariance allows for the calculation of the aver- 
ages lDll without a direct integration |3^| . It can be shown 



341 that the tensor Ti. 



must be proportional to a 



sum of products of Kronecker's deltas symmetric with 
respect to the indexes of the same type (i.e. the i's or 
the j's); the normalisation coefficient can then be fixed 
by comparison with the case R^ = I Vfc. For r = 1 



therefore Z(R) = 



tr(R) 
N ' 



(D2) 



For the case r = 2 one similarly obtains 
Tijki oc (6ij6ki + SuSjk)- 
Thus, the integral IDll can be written using only traces 
tr(i?( 1 ')tr(i?( 2 ))+tr(i?( 1 )i?( 2 )) 



i? (i) i? (2) T 



N 2 + N 



(D3) 



If the operators are simultaneously diagonalisable, it 
is easy to see that the integral depends only on their 
spectrum. If R is a generic unitary operator, setting 
I(R 1 RJ) — 12{R) and A = one obtains that 



1 2 (R) = 1-A 



N 2 



(D4) 



When R has only two eigenvalues, e zXl and e tX2 , with 
the same multiplicity N/2 (i.e., when it is a generic lQ 
operator), the previous formula reduces to 



('o 1 /x 2 )= 1 "^ sin2 (^ 1 )- 



(D5) 



It is interesting to calculate the same integral in a dif- 
ferent way. By introducing the probability < p < 1 of 
being in the first eigenspace, the state can be written as 



l^> = vWO + v^-pIV*). 

where the normalised state | ipj ) is an eigenstate with 
eigenvalue e lXj , and all the unnecessary phases have been 
reabsorbed. The integral in formula |D5I becomes then 



( e lA i 
V e iX '< 



pe lM + (1 -p)e tA2 \ V{p)dp 



= i _ 4s in 2 (^) • [<p)(l - (p)) ~ a 2 p ] . (D6) 

A freedom in the choice of the probability density 
function V(p) corresponds to a freedom for the measure 
dfi{ip). Due to the symmetry of Ai and A2 in the uniform 
distribution, (p) is necessarily 1/2. Formulae ID5l and lD6l 
then coincide when 



, I- A 



1 



4(1 + N) 



In the framework of quantum computation, the dimen- 
sion of the space is exponential in the Q.C. size, N = 2 rlq , 
so cr 2 is exponentially small. This behaviour, known in 
literature as the concentration of measure phenomenon 
[2ll |. is much more general than the case studied here: 
in a space with a large number of dimensions, the de- 
viation of a Lipschitz function from its average value is 
extremely small. Without entering into the details, in 
general V(\I — (I) \ > t) ~ e~ aNt (with a a constant), 
so that the deviation is of order 0(l/>/~N). In other 
words, almost all the vectors are "typical" with respect 
to Lipschitz functions with the uniform measure. 

The generalisation of formula ID5I for an operator with 
many different eigenvalues is quite complicated. The 
limit of small errors has however a simple interpretation 
in terms of the variance of the eigenvalue phases. For 
generic eigenvalues {e lXl , . . . , e lAjv } of R, let 



(A) 



N 



and 



N 



(D7) 



The small errors limit corresponds to eigenvalue phases 
with small spread (ax <§C 1); in this limit the trace |tr(i?)| 
can be approximated by retaining only the terms up to 
the second order, 



iX ; 



i[\,-(A>] 



1 



7V 2 (1- 



so that formula ID4I gives a value of the integral which 
depends only on the variance of the eigenvalue phases: 



MR) = 1 



Aa\ 



0(a 



(D8) 
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Problem EU can be generalised to involve more vec- 
tors with correlated distributions, for instance a variable 
| ip ) with uniform distribution dfi(ip) in a generic vector 
space, and another variable \<j>) with uniform distribu- 
tion d/j,(<p\^) in the subspace orthogonal to \ ip), like in 



J 



R (r) 



= / d^wm) n(V'ii? (fc) i0>: 



where r is a number of operators and the symbol * W 
means either a null or a conjugation sign (i.e. some her- 
mitian products can be conjugated). In reference [34| it 
is shown that the number of conjugations must be equal 
to that of "non conjugations" in order to have a non-zero 
value for the integral. The J integrals can be written as 
functions of the X integrals; as an example, the case r = 2 
will be worked out here. The first step consists in insert- 
ing a unitary operator S<p mapping a fixed vector | ) 
into \ip). Changing the integration variable from | <f> ) to 

= | (f> ) allows to express the measure in a simple 

way, because d/i (</>') becomes simply the uniform measure 
in the (N — l)-dimensional subspace Span {| i ) | i ^ 0}: 

MR) = j( R , J=/ d^)dv{<t>\ip)\(ip\R\4>)\ 2 

d^){^\RS^[ J <t>){<t> \Si]sl&\il>) 

dK1>) ( i> \R% [ / dfi(<t>') 1 $ ) < <t>' |] %W | v > • 



It follows from result lD2l that f dfJ,((p')\ ft ){ (j>' | can 
be replaced by -^j-Ijv-i, where Ijy-i is the identity in 
the appropriate subspace with N — 1 dimensions. The 
action of maps this identity to the projector onto 
the subspace orthogonal to \ip), that is S^In-iSI, = 
I — | tjj ) { ij) | , therefore 



MR) = 



1 



N- 1 
1 

TV — 1 



dni^i^lRS^lN^Sl&^j) 

d / i(V)(^|i?fi-|^)(vili? t l^). 



Splitting the integral and inserting the result ID8I for 
T2{R) shows once again that, in the limit of unitary R 
close to the identity, the result depends only on the di- 
mension of the space and the spread of the eigenvalue 
phases o\ : 



MR) 



1{RR))-1 2 {R) _ N 
N - 1 ~ N 2 - 



+ 0(4). (D9) 
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